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We considered the interaction of semiconductor quantum register with noisy environment leading 
to various types of qubit errors. We analysed both phase and amplitude decays during the process of 
electron-phonon interaction. The performance of quantum error correction codes (QECC) which will 
be inevitably used in full scale quantum information processors was studied in realistic conditions in 
semiconductor nanostructures. As a hardware basis for quantum bit we chose the quantum spatial 
states of single electron in semiconductor coupled double quantum dot system. The modified 5- 
and 9-qubit quantum error correction (QEC) algorithms by Shor and DiVincenzo without error 
syndrome extraction were applied to quantum register. 5-qubit error correction procedures were 
implemented for Si charge double dot qubits in the presence of acoustic phonon environment. \~ 
matrix, Choi-Jamiolkowski state and measure of decoherence techniques were used to quantify 
qubit fault-tolerance. Our results showed that the introduction of above quantum error correction 
techniques at small phonon noise levels provided quadratic improvement of output error rates. The 
efficiency of 5-qubits quantum error correction algorithm in semiconductor quantum information 
processors was demonstrated. 



I. INTRODUCTION 



In recent years much attention is attracted to the influence of decoherence on quantum communication channels [TJ- 
[3] and quantum information processing [3 6 J . Measure of decoherence approach [7J |Hj occurred to be helpful tool for 
quantitative evaluation of quantum state distortion due to noisy environment, x-matrix and Choi-Jamiolkowski state 
representations are very efficient for quantum channel description and study PH31 15HT2] . In this paper we combine 
these approaches to analyse the influence of error correction |13H19j on quantum registers. The representative example 
of solid-state qubit - electron in semiconductor double dot interacting with phonons is also considered. 



II. DESCRIBING QUANTUM CHANNELS 



We consider a qubit system described by a density matrix p of size d x d. By definition, the operator p is Hermitian, 
positive scmidcfinite and has trace one QJI31 [23 121] ■ An evolution of a density matrix is described by a quantum 
operation £ (also called a quantum channel or a stochastic map) PJEJ 1514X2] ■ We use these terms interchangeably. 
As usual a linear map £ is called a quantum operation when it preserves trace (in Schrodinger picture) and is 
completely positive [Tj. In our investigation, a qubit system is placed in an environment that acts on these qubits 
independently. That means that the quantum operation describing the impact of the environment on the system of 
qubits is separable. Choi has shown 22 that every quantum channel can be represented as a sum of Kraus operators 
(Kraus representation [23]): 

n n 

s[p]=Y,EipEl J2^A=i- a) 

i=l i=l 

This statement is also a sufficient condition for a map to be a quantum channel. Equation shows a way to describe 
a channel. 

There are several ways to describe a quantum channel mathematically. In the sake of convenience, we will use a 
Choi-Jamiolkowski state |24) r in simulation of quantum circuits and a ^-matrix representation in the calculation of 
measure of decoherence. Below we consider both of these representations. 
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A. x _ma trix representation 

Every linear map £ which is a quantum operation can be written in the form of: 

d 2 -i 



£[p}= Y X»pE a pE i 

a, 0=0 



(2) 



where the set of {E a }^_Q matrices is the basis in Wf 2 space. 

The coefficients Xap form \~ matrix of dimension d 2 x d 2 3J. From the construction of \ it is seen that the matrix 
is Hermitian. Since = £[p], we have 



d 2 -l \ rf 2 -l 

K a,P=0 J a,p=0 



d 2 -l 



Y X* a p E ppE\ = Y XapEapEp = Y Xp a Epp& a . Thus X p a = X* a 



(3) 



a, 0=0 



a,P=0 



Another limitation is imposed by the fact that channel £ preserves trace because it is a quantum operation. Since 



Tr [Ep] = Tr T,t,pt X a pE a pE\ = Ea )/3 =o Tr 



X a ppE^E a 



= Tr 



P J2 a p=o XapEiE a = Ti[p\ , we obtain 



d*-i 



Y X a pE\E a = I. 



(4) 



a, 0=0 



Since %-matrix of 4 x 4 dimension is Hermitian, it can be defined by 16 real numbers. The limitations due to 
trace conservation shown in Eq. Q allows us to describe \ with 12 real parameters, ^-matrix formalism provides a 
simple method of verification of complete positivity. A map £ is completely positive if and only if ^-matrix is positive 
semidefinite [231. 



B. Choi-Jamiolkowski state r representation 



Operator x i s given by 



-l 



X=Y X a p\E a ))((Ep\, (5) 

a, 0=0 

where \A)) = [An, A 1( i, A 2 i, A 2d , A d i, A dd ] T is a supervector in Liouville space (L-space). Supervectors 
{l^a))}^ 1 form the basis in n d 2 - 



Choi-Jamiolkowski state is defined as 



= (e®i)\a d2 )(si d2 \, (6) 



where the state jO^a) = £i=o i s a maximally entangled state. It is worth noting that ^1-0) = \^d 2 ) an d 

\A)) = (A<g)fj \I)). These properties are followed from the above definitions. The quantum operation is uniquely 
determined by the Choi-Jamiolkowski state. This fact is known as Choi-Jamiolkowski isomorphism. Let us show that 
t = x/d: 

d 2 -l . 

x= Y x*e(Ea®f)\i)){{i\(Ep®i) =(s®fyi)){{i\ = d(8®i)\n d2 ){n d 2\ = Td. (7) 

a,p=0 



3 



III. MEASURE OF DECOHERENCE 

To estimate the quality of quantum correction we use the concept of measure of decoherence D [6HH1 EH) • By the 
definition, measure of decoherence is the maximum over all the states of the operator norm of matrix p out — p m : 



D 



SUp ||p ou t — Pin 1 1 , 

Pin 



(8) 



l a l PQ- By spec(A) denote the 



where the operator norm of Hermitian matrix A is given by ||^4|| = max nespcc ( j 4) 
spectrum of operator A. 

Let us consider a qubit (d — 2). In this article it is shown that if the influence of the environment has the Kraus 
representation in Pauli basis with ^-matrix in the form of 



X 



/Xo 

Xi o 

X 2 o 

Vo o o X3/ 



then the expression for measure of decoherence is simplified 

D = max{xi +X2,Xi +X3,X2 + X3> = Xi + X2 



X3 



l{Xl,X2,X3}- 



(9) 



(10) 



Let Xmk = Xlnk + *Xmfc> then usin S Eq. [3] we get the relations x\l k = Xfm^Xmk 
system of 4 equations: 



-xll 



Using Eq. [4] we obtain a 



Xoo 


= 1 - 


im 
Xl2 


= X30 


v im 


= X2 C 


im 
^ A23 


= X r i c o 



X22 — X33 



Consider an arbitrary ^-matrix of 4 x 4 dimension: 



A - Xi ! - X r 2 c 2 . - X r 3 e 3 X r o e i + ixS? Xo! + iX$ Xhl + *X m \ 



X01 - 1X01 
Xo! - iX$ 
X55 ~ 



xii Xia + *Xo3 xi 3 -*Xo2 

X r i2 - iXol X22 X23 + «Xo! 

XS+ftS X^Xoi X?3 J 



Let us rewrite it in other variables 



- Xl - X2 - X3 X4 + «X5 X6 + iX7 X8 + «X9 ' 

X4 - iX5 Xi X10 + «X8 X11 - «X6 

X6 - iX7 X10 - «X8 X2 X12 + iX4 

X8 - «X9 Xll + «X6 Xl2 - «X4 X3 

In general, an arbitrary density matrix can be written as p = + £P ^\ 
herence depends on the components of ^-matrix as follows: 



(11) 



(12) 



(13) 



2. Then measure of deco- 



D( X ) = max 2 X4 P X + 2 Xe P y + 2 X %P Z ± 

± [(XlO + Xll + X2 + X3 + X7 + xl + 2X2X3 - 2xiiX7 + 2xioX9) Pl + 
+ (X10 + XI2 +xl+xi + xi+xl + 2xiX3 - 2xioX9 + 2xi2X 5 ) P*+ 
+ (Xn + X12 + xl + xl + xl + Xi + 2xiX2 - 2xi2X5 + 2 X u X t) P*+ 
+ 2 (2xioX3 + X10X2 + X12X7 + X5X7 + X10X1 + X9X1 - X11X12 - X11X5 - X2X9) P x P y + 

+ 2 (2Xl2Xl + X12X2 + X12X3 + X2X5 + XllX9 + X7X9 - XloXll - X10X7 - X3X5) Py P z + 

1/2 

+2 (2xnX2 + X11X3 + X10X5 + XlXll + X3X7 + X5X9 - X10X12 - X12X9 - XlX7) PxPz] 



(14) 
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FIG. 1: An n-qubit QECC and the scheme for calculating r-matrix in case of arbitrary error correction, n-qubit QECC is 
shown on this Figure. \Q) is a maximum entangled two-qubit pure state. The first qubit of the state \Q) remains unchanged 
during the process of QEC, the second qubit of \Q) is the first qubit of n-qubit QECC. 



\Q) 



Noise 



FIG. 2: Scheme for calculating measure of decoherence in case of an error in a qubit. The first qubit remains unchanged, the 
second is affected by a quantum channel. 



In the case Xi — Xe — Xs — 0, the square of the measure of decoherence is a quadratic form, which can be reduced 
to canonical form by introducing the new variables Q x , Q y , Q z : 

D 2 = max [axQl + a 2 Q 2 y + a 3 Q 2 z ] . (15) 

Let us consider the function / = a\Q 2 + a 2 Q 2 + CL3Q 2 , where the coefficients {ai, a 2 , 03} arc the eigenvalues of the 
quadratic form /. With this change of variables: 

Ql + Q 2 y + Ql<i 
Ql>0,Q 2 y >0 7 Q 2 z >0 (16) 

/>o. 

Since for all Q x ,Q y ,Qz the function / should be nonnegative, we get ai > 0, a 2 > 0,a 3 > 0. Hence the function 
/ is nondecreasing, the maximum of this function is achieved on the boundary of Eq. (16 1. The point {Q x ,Qy,Qz) 



of space at which there is the maximum of measure of decoherence is limited to the set{(l,0,0), (0, 1,0), (0,0, 1)}. 
Therefore, we have 

D =v/max{spec(/)}. (17) 



IV. SIMULATION OF NOISY QUANTUM CHANNELS AND QEC 

To find an analytical expression for the measure of decoherence it is necessary to obtain x-matrix, that describes the 
influence of an environment on the qubit. The given \- matrix depends on parameter q, it is the probability of the qubit 



to decohere due to interaction with an environment. The function D{p) can be calculated using Eqs. ( 10 1 , ( 14 ) and ( 17 1 . 
Despite the usability of xip)i the calculation of this function according to definition of ^-matrix j n complicated. It 
is difficult to write a quantum operation in the form of Eq. |2| in case of simulating QECC. This problem can be 
avoided by calculating Choi-Jamiolkowski state r. The scheme for calculating r-matrix in the case of QEC is shown 
on the Fig. [l] It is easy to find \ using r. From Eq. (5| we can derive Xap coefficients which form ^-matrix. 
The above procedure can be summarized as the following sequence: 

1. We write down the evolution of the density matrix for a single qubit interacted with an environment as an 
operator-sum. This decomposition depends on the parameter q. The parameter characterizes the amount of 
errors. 

2. We calculate r(q) state after the influence of an environment (Fig. [2]) . We obtain x{l) state using r(q). 



3. From x(q) we § e t x(9) _ma t r i x - 
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FIG. 3: Scheme for calculating measure of decoherence in case of bit flip correction. 



4. We calculate measure of decoherence D (q) for qubit in an environment. Here D (q) is a probability for p to 
decohere. 

5. We replace q with p in the operator-sum based on the equality Do(q) = p. 

6. We obtain the limitation on p based on the positivity of x(p)-matrix. 

7. Next we consider the error correction code. We take N different values of p, using the existing limitations on 
the probability p. We calculate Choi-Jamiolkowski state r(p) after applying QECC. We get x(p) state using 
T(p). 

8. From x(p) we get x(p)-matrix. 

9. We calculate measure of decoherence D(p) of qubit after interaction with an environment and quantum correc- 
tion. For n different p we take n values of D, n is the number of qubits in QECC. 

10. Since the quantum correction code gives a polynomial in p improvement we have 

n 

D(p)=J2<*iP i - (18) 

We solve n equations and calculate coefficients {ai}f =1 of the polynomial of the n-th degree. Hence we establish 
polynomial expression for D. 

Next we consider the results of the algorithm for calculating D(p). 



A. Bit and Phase Flips 

Let bit flip environment changes |0) to |1), and |1) to |0) with probability p' . This action can be conveniently 
written as Eq. (JiJ) with the Pauli matrices {I , X ,Y , Z} as the basis {E^I^Jq 1 for d = 2 (single qubit). Then the 
operator-sum is written as: 

£[ P ] = {l-p')ipi + p l XpX. (19) 

This representation in the form of Kraus decomposition allows us to write %-matrix of the quantum channel as a 
diagonal matrix in Eq. ^ with parameters xo — 1 — p' ■> Xi — p' an d X3 = X4 = 0. It follows from Eq. (10) that 
measure of decoherence D = p = p'. Then ^-matrix f the quantum channel in case of bit flip is equal to 

'1—p ON 



pOO 

0/ 



Xo=| J 00 |. (20) 



According to %-matrix positivity criterion, bit flip channel is completely positive if and only if < p < 1. 

Let us consider a bit error correction code based on a decoding of type called majority voting Scheme for finding 
measure of decoherence in the case of majority voting is represented in Fig. [3] 

^-matrix is calculated following the above algorithm: 




ON 

3p 2 - 2p 3 



0/ 



(21) 
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FIG. 4: Scheme for calculating measure of decoherence in case of phase flip correction. 
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FIG. 5: 5-qubit DiVincenzo-Shor QECC. R is the unitary operator [18, 27 . The code does not contain any measurement 
operations. 



According to Eq. (10), the measure of decoherence 

D(p) = 3p 2 - 2p 3 . 



(22) 



The procedure of finding the result of an action and the result of a correction in case of phase flip is similar to the 
case of bit error correction. The only difference is the modified scheme of majority voting, shown in Fig. [4] 

B. Depolarizing Channel 

It is convenient to write an operator sum in case of the depolarizing channel in the form of Kraus decomposition: 



i 



-p) Ipl + 1 (x P x + YpY + ZpZ] 



2 V 



(23) 



Since all the matrices in the decomposition are the Pauli matrices, the matrix \o is diagonal and is equal to 

(l- \p Q \ 



Xo 



p/2 
Op/2 
\ p/2/ 



(24) 



The parameters in the decomposition in Eq. 23 are selected to satisfy Dq — p. From positivity of xo _ma trixp € [0, 2/3]. 

Let us consider the result of the 5-qubit DiVincenzo-Shor QEC [TU [TB] , shown in Fig. |5j after the action of the 
depolarizing channel. In the case of the depolarizing channel it is possible to obtain a general view of x-matrix, which 
describes the action of the error correction code: 



X = 



3xi 














Xi 














Xi 














Xi 



(25) 



where 



Xi = p 2 (15 - 50p + 60p 2 - 24p 3 ) / 2. 



(26) 
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FIG. 6: 9-qubit Shor QECC. This code is a composition of bit and phase flip error correction codes shown in Figs. [3] and [§ 
The code does not contain any measurement operations. 
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FIG. 7: QEC in case of a depolarizing channel. Measure of decoherence D vs. the probability of an error in a qubit p. The 
results are consistent with the numerical simulation 1271. 



It follows from Eq. ( 10 1 that measure of decoherence is equal to 

Dip) = P 2 (15 - 50p + 60p 2 - 24p 3 ) 



(27) 



Let p — »• 0. For 9-qubit Shor code [T3] shown in Fig. [6] measure of decoherence has the form 

Dip) = 36p 2 . (28) 
Combining Eqs. (27) and (28 1, we obtain the result of QEC shown in Fig. [7] in case of the depolarizing channel. 



Let us consider the basis 
is written in the form of 



V. QEC IN SI DOUBLE DOT CHARGE QUBITS 
A. Amplitude Damping 

1+) = ^ [|0) — |1)]. In this basis the evolution of the density matrix 



l-p__(0)e- rt /9+ _(0)e- r '/ 2N 
p_+(0)e- rt / 2 p__(0)e- rt 1 



(29) 



denote by T(t) a relaxation rate. Let us rewrite Eq. (|29[) to obtain the representation in the form of an operator-sum: 



£[p} = 



i 

e 







1 







rt/2 P[ Q e -rt/2 I + 



VI - e" 



r/ 











P 



VI - e" 



17 











(30) 



Let us consider the state r. The nonzero eigenvalues of the state are (1 ± e _r *)/2. Since Ft > 0, we see that 
all eigenvalues of the matrix r are nonnegative, that is r is positive semidefinite. Therefore, the map in Eq. (30) is 
completely positive. 

We represent the quantum operation shown in Eq. (30) in the form of xo _ma trix in Pauli matrix basis: 

, 2 



,-rt 



Xo 



) \ 



(e~ rt -l) i(e- rt -l) 

-i(e- rt -l) (e- rt -l) 

V (l~e- rt ) (l-e- rt / 2 )7 



(31) 



According to Eq. ([8j) measure of decoherence is equal to D n = 1 — e 
that measure of decoherence in case of the 5-qubit QECC 

D(p) = 5p 2 (3-3p + p 2 )/8. 



p. It follows from the proposed algorithm 



(32) 



It can be concluded that for small p the efficiency of the error correction in the case of amplitude damping is 8 times 
higher than in the case of the depolarizing channel. 



B. Phase Damping 



The evolution of the density matrix in a single operation can be written as 



Poo(0) Poi(0)e- s ' : 
p 10 (0)e- B2 Pll (0) 



denote by B (t) a spectral function. Using representation in Eq. (33), we can write the operator-sum: 



(33) 



E\p] = 



-B 2 /2 



-B 2 /2 



e- B2 / 2 
e~ B ' 2 ' 2 



A/1 - e- B ' 2 0^ fVl-e- 







o VT 



o VT 



,-B< 



(34) 



Consider Choi-Jamiolkowski state r. The nonzero eigenvalues of the state are ^1 ± e B ^ / 2. Since B 2 > 0, we 

have nonnegative eigenvalues of matrix r, that is the state r is positive semidefinite. Hence, the map in Eq. (34) is 
completely positive. 

Let us represent the quantum operation of Eq. (1341) in the form of xo -ma trix in Pauli matrix basis: 



Xo 



,-B 






! ) /2 






\ 



1 



(35) 



(l-e^)/2; 
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FIG. 8: QEC in Si double dot charge qubit. Measure of decoherence D vs. the computational cycle time t. The figure shows 
the error correction after N=68 operations. The result are consistent with the numerical simulation [27| . 



It follows from Eq. ( 10 1 that measure of decoherence is equal to Dq = yl — e B 
algorithm we calculate measure of decoherence in case of the 5-qubit QECC 

D(p) = I0p 2 (l-2p + p 2 ). 



2 = p. According to the proposed 



(36) 



It can be concluded that for small p the efficiency of the error correction in the case of amplitude damping is 1.5 times 
higher than in the case of the depolarizing channel. 



C. QEC Results 

During calculations, transmission or information storage there occurs amplitude and phase damping. To fight the 
process of decoherence one can use error correction algorithms, which have already been discussed in this paper. 

Let us consider an error correction in silicon, where the qubit represents as a state of an electron in the double 
quantum dot 5 , 28 - 30 . To describe the entire system of five qubit we use our results for measure of decoherence in 
Eqs. (32 1 and (36). To determine relaxation rate and spectral function, we use the following formulas 



E 2 k 3 
4irps 2 h 



exp(-a 2 fc72) 1 



sin(/cL) \ 
kL J ' 



(37) 



B 2 (t) 



~2 



TT 2 hpS 3 



2 f . sin 2 ( g Lcose)exp(-aV/ 2 ) . 2 qst 
q da j smBdB sin , 

J q 2 

o 



(38) 



where deformation potential 5 = 3.3 eV, speed of sound s = 9.0- 10 3 m/s, crystal density p = 2.33 g/sm 3 , the distance 
between points L = 50 nm, radius of points a = 3 nm. 

Depending upon the quantum gate operation required at current algorithmic step, a qubit experience phase or 
amplitude damping error during quantum computing on double dot qubit. Hence we define measure of decoherence 
in this case as the maximum of the two measures. Combining Eqs. (32 1 and (36), we obtain 

(39) 
(40) 



D (t) = max{pi,p 2 }) 

D(t) = max{ 5^(3 - 3 Pl +p\)/% , 10^(1 - 2p 2 +p%) } , 

where pi = 1 - e r( - t)t and p 2 = (l - er B{t)2 ^j j 2. 

Let quantum computer perform N operations. By p denote the probability of an error during one cycle time. 
Substituting p\ for Np\ and p 2 for Np 2 in Eqs. (39) and (40), we obtain the result of QEC shown in Fig. [8] 
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VI. CONCLUSIONS 

By using x-matrix, Choi-Jamiolkowski state, and measure of decoherence techniques we analysed the influence of 
noise on quantum bits. The analytical expressions of measure of decoherence were obtained for practically important 
wide subset of quantum channels. We have shown that the introduction of DiVinccnzo-Shor quantum error correction 
algorithm is helpful for double dot charge qubits and reduces effective error rate drastically. 
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